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Abstract 

Compressible isothermal turbulence is analyzed under the assumption of homogeneity and in 
the asymptotic limit of a high Reynolds number. An exact relation is derived for some two-point 
correlation functions which reveals a fundamental difference with the incompressible case. The 
main difference resides in the presence of a new type of term which acts on the inertial range 
similarly as a source or a sink for the mean energy transfer rate. When isotropy is assumed, 
compressible turbulence may be described by the relation, — |£ e gr = J- r (r), where T r is the radial 
component of the two-point correlation functions and e c g is an effective mean total energy injection 
rate. By dimensional arguments we predict that a spectrum in /c~ 5 / 3 may still be preserved at small 
scales if the density-weighted fluid velocity, / o 1 / 3 u, is used. 

PACS numbers: 47.27.eb, 47.27.ek, 47.27.Gs, 47.40.-x 
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Introduction. Fully developed turbulence is often seen as the last great unsolved problem 
in classical physics which has evaded physical understanding for many decades. Although 
significant advances have been made in the regime of wave turbulence for which a systematic 
analysis is possible the regime of strong turbulence - the subject of this Letter - continues 
to resist modern efforts at solution; for that reason any exact result is of great importance. 
In his third 1941 turbulence paper Kolmogorov derived an exact relation for incompressible 
isotropic hydrodynamics in terms of third-order longitudinal structure function and in the 
asymptotic limit of a high Reynolds number (-Re) |2j. Because of the rarity of such results, 
;he Kolmogorov's universal four-fifths law has a cornerstone role in the analysis of turbulence 
3)]. Few extensions of such results to other fluids have been made; it concerns for example 
scalar passively advected such as the temperature or a pollutant in the atmosphere, quasi- 
geostrophic flows or astrophysical magnetized fluids described in the framework of (Hall) 
MHD j4J. It is only recently that an attempt to generalize such laws to axisymmetric 
turbulence has been made but an additional assumption is made about the foliation of the 
correlation space 

The previous results are found for incompressible fluids and to our knowledge no universal 
aw has been derived for compressible turbulence (except for the wave turbulence regime 
6j) which is far more difficult to analyze. The lack of knowledge is such that even basic 
statements about turbulence like the presence of a cascade, an inertial range and constant 
flux energy spectra are not well documented 0]. That is in contrast with the domain 
of application of compressible turbulence which ranges from aeronautical engineering to 
astrophysics js-lQ|- In the latter case, it is believed that highly compressible turbulence 



llf whereas in the former case Re is relatively 



controls star formation in interstellar clouds 
smaller. 

In that context, the pressure- less hydrodynamics is an interesting model to investigate 
the limit of high Mach number compressible turbulence whose simplest form is the one 
dimension Burgers equation which has been the subject of many investigations 12j. Among 
the large number of results, we may note that with exact field-theoretical methods it is 



possible to find explicit forms of some probability distributions 



13j; it is also possib 



derive the corresponding exact Kolmogorov law for the third-order structure function 



e to 
3|- 



In the general case, our knowledge of compressible hydrodynamic turbulence is mainly 
limited to direct numerical simulations jl^]. The most recent results for supersonic isother- 
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mal turbulence with a grid resolution up to 2048 3 [15| reveal that the inertial range velocity 
scaling deviates substantially from the incompressible Kolmogorov spectrum with a slope 
of the velocity power spectrum close to —2 and an exponent of the third-order velocity 
structure function of about 1.3. Surprisingly, the incompressible predictions are shown to be 
restored if the density-weighted fluid velocity, p 1 / 3 u, is used instead of simply the velocity u. 
Although a —2 spectrum may be associated with shocks - like in one dimension - it seems 
that their contribution in three dimensions (3D) is more subtle. Generally speaking it is 
fundamental to establish the equivalent of the 4/5s law for compressible turbulence before 
going to the more difficult problem of intermittency 

In this Letter, compressible isothermal hydrodynamic turbulence is analyzed in the limit 
of high Re. We shall investigate the nature of such a compressible turbulence through an 
analysis in the physical space in terms of two-point correlation functions. In particular the 
discussion is focused on the isotropic case for which a simple exact relation emerges. The 
theoretical predictions illuminate some recent high-resolution direct numerical simulations 
made in the astrophysical context. 

Homogeneous compressible turbulence. We start our analysis with the following 3D com- 



pressible equations 



d t p + V- (pu) = 



d t (pu) + V • (puu) = -VP + fj.Au + ^V(V ■ u) + f , 

where p is the density, u the velocity, P the pressure, p the coefficient of viscosity and f a 
stationary homogeneous external force acting at large scales. The system is closed with the 
isothermal equation P = C 2 s p where C s is the speed of sound. The energy equation takes 
the form 

d t (E) = -p((V x u) 2 ) - ^((V • u) 2 ) + F , (1) 

with ( ) an ensemble average (which is equivalent to a spatial average in homogeneous turbu- 
lence), E = pu 2 /2 + pe the total energy, e = C 2 ln(p/p ) (po is a constant density introduced 
for dimensional reasons) and F the energy injected. 

The relevant two-point correlation functions associated with the total energy may be 
obtained by noting that for homogeneous turbulence 

(5(pu)-6u) = 2<pu 2 )-<(p + p>.u'>, (2) 

(5p5e) = 2(pe)-(pe' + p'e), (3) 



where for any variable £, 5C, = £(x + r) — £(x) = £' — £. Then, we find 

^ + 2 n ~ r) = (E) - \ (5(pu) • 6u) - \ (5p5e) , (4) 

where K(r) = (pu ■ W/2 + pe') = (R) and K(-r) = (p'u' ■ u/2 + p'e) = (R). Note that for 
homogeneous compressible turbulence the relation, lZ(r) = TZ(—r), holds only when isotropy 
is assumed whereas it is always valid in the incompressible limit for which 1Z is reduced to 
a second-order velocity correlation function [l8|. As we will see below, relation (|4|) is very 
helpful for deriving an exact relation for some two-point correlation functions. In practice, 
we shall derive a dynamical equation for 7£(r); first, we have to compute 

d t (pu ■ u') = (pu • d t u' + u' ■ d t (pu)) 

= (pu ■ (-u' ■ W - ivP0> (5) 
+{a! ■ (-V ■ (puu) - VP)} + 2V + , 

where for simplicity TD and 2J 7 denote respectively the contributions to the correlation of 
the viscous and forcing terms. By remarking that 

(^u ■ VP') = (pC 2 s u^) = (pu e d' e e') = (V • (pe'u)) , 

and that 

(u' • V'(pu • u')> = (V • (p(u • uV) - p(u • u')(V • u')) , 
we can rewrite ([5]) in the following way 

d t (pu-u') = (-W ■ V'(pu ■ u') - V • (pe'u)) 

-(V • (p(u ■ u')u + Pu')) + 2V + 2F 
= V r • (-p(u ■ u')5u + Pu' - pe'u) 

+(p(u-u')(V' -u')) + 2V + 2T. (6) 

Secondly, we have to complete the computation with 

d t (pe') = (pcV + e'd t p) = (C s 2 4 d t p' + e'd t p) 

P 

= (-C s 2 4V-(p'u')-e'V-(pu)} 
P 

= -(V'-(C s 2 pu')) + (p'u'-V'(C 2 -^)) 

-(V-(pe'u)}. (7) 



By noting that 



(P'u'ed'^)) = -(V ■ (pe'u')> + (e'V • (pu')> 



we obtain after simplification 

dt(pe') = V r • {-pe'Su - Pu') + (pe'(V ■ u')) . (8) 
The combination of and (jSJ) leads to 

3tft(r) = ((V' • u')R) +V + J r (9) 
+V r • (^—RSu — 2^ >u ' ~~ 2 /0e ' u 

The same type of analysis may be performed for TZ(—r) which eventually leads to the 
dynamical equation 



±<(V ■ u')i?) + i((V • u)^} + + p + ^ + 

+^V r ■ (-{R + R)Su - i(Pu 7 - pe'u - P'u + peu 

where X> and J 7 denote respectively the additional contribution of the viscous and forcing 
terms. 

Local turbulence. For the final step of the derivation we shall introduce the usual as- 
suction SP ec lflc t o 3D My deve.oped turbul ence with a di.ec. e„e rgy c as cade Q. in 
particular, we suppose the existence of a statistical steady state in the infinite Reynolds 
number limit with a balance between forcing and dissipation. We recall that the dissipation 
is a sink for the total energy and acts mainly at the smallest scales of the system. Then, 
far in the inertial range we may neglect the contributions of T> and T> in equation ([10 
The introduction of structure functions leads to the final form 

- 2e = ((V • u')(R - E)) + ((V • u){R - E')) (11) 



19]. 



+ V r - 



5(pu) ■ 5u 

- + 5p5e - C;5p 



Su + SeS(pu 
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where SX = (X + X')/2 and e is the mean total energy injection rate (which is equal to 
the mean total energy dissipation rate; see relation (CQ)). Note that at relatively small r the 
function, R — E (and R — E'), and its derivative are negative since the correlation between 
two points is maximum if the points are the same. 



It is straightforward to show that in the limit of incompressible turbulence we recover 
the well-known expression; indeed we obtain (with p — > p Q — 1) 

_4e = V r • ((5u) 2 5u) , (12) 

which is the primitive form of the Kolmogorov's law. An integration over a ball of radius r 



leads to the well-known expression 20|, — (4/3)er = ((6u) 2 5u r ) , where r means the radial 
(often called longitudinal) component, i.e. the one along the direction r. 

Isotropic turbulence. Expression fllip is the main result of the Letter. It is an exact re- 
lation for some two-point correlation functions when fully developed turbulence is assumed. 
It is valid for homogeneous - non necessarily isotropic - 3D compressible isothermal turbu- 
lence. Note that the pressure contribution appears through the term C 2 6p and is therefore 
negligible in the large Mach number limit (C s — > 0). When isotropy is additionally assumed 
this relation can be written symbolically as 

- 2e = S(r) + ^d r (r 2 T r ) , (13) 

where T r is the radial component of the isotropic energy flux vector. In comparison with 
the incompressible case ffl2|) . expression (TIB"]) reveals the presence of a new type of term S 
which is by nature compressible since it is proportional to the dilatation (i.e. the divergence 
of the velocity). This term has a major impact on the nature of compressible turbulence 
since as we will see it acts like a source or a sink for the mean energy transfer rate. Note 
that S consists of two terms which account for two-point measurement approach. 

Discussion. We may further reduce equation ffl3|) by performing an integration over a 
ball of radius r. After simplification we find the exact relation 

2 1 r 

-sr = — S(r)r 2 dr + F r (r) . (14) 

3 r z J 

We start the discussion by looking at the small scale limit of the previous relation which 
means that the scales are assumed to be small enough to perform a Taylor expansion but 
not too small to be still in the inertial range. We obtain S(r) = S(0) + rd r S(0) = rd r S(0) 
which leads to 



£ + lrd r S{0) 



2 

--Sesr = J>(r) . (15) 
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Note that we do not assume the cancellation of the first derivative of S at r = although 
the function, R — E, reaches an extremum at r = 0; the reason is that this function is 
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weighted by the dilatation function which may have a non trivial form. We see that at the 
leading order the main contribution of S(r) is to modify e for giving an effective mean total 
energy injection rate £ e g. Then, the physical interpretation of (j!5p is the following. When 
the flow is mainly in a phase of dilatation (positive velocity divergence), the additional term 
is negative and e e g is smaller than e. On the contrary, in a phase of compression d r S(0) is 
positive and e e ff is larger than e. 




FIG. 1: Dilatation (left) and compression (right) phases in space correlation for isotropic turbu- 
lence. In a direct cascade scenario the flux vectors T (dashed arrows) are oriented towards the 
center of the sphere. Dilatation and compression (solid arrows) are additional effects which act 
respectively in the opposite or in the same direction as the flux vectors. 

An illustration of dilation and compression effects in the space correlation is given in 
Fig. [H In both cases, the flux vector T (dashed arrows) is oriented towards the center 
of the sphere (r = 0) since a direct cascade is expected. Dilatation and compression act 
additionally (solid arrows): in the first case, the effect is similar to a decrease of the local 
mean total energy transfer rate whereas in the second case it is similar to an increase of the 
local mean total energy transfer rate. 

The discussion may be extended to the entire inertial range (i.e. for larger values of r) 
when the (turbulent) Mach number is relatively high. In this case the analysis is focused on 
expression (|T4|) for which we have already noted that a term like R — E is mainly negative. 
It is interesting to note that S(r) is composed of two types of term which are different by 
nature. First, there is the dilatation dominated by the smallest scales in the flow - the 
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shocklets - which mainly give a negative contribution with a fast variation 22[. Secondly, 
there is the correlation R — E which derives most of its contribution from relatively larger 
scales with a slower variation. This remark may lead to the assumption that both terms are 
relatively decorrelated 19||. Then S(r) may be simplified as (by using relation fll])) 

(16) 



S(r) ~ -^J(V-u) 



-5(pu) ■ 5u + -Sp5e 



The previous expression is not derived rigorously but it may give us some intuition about 
its contribution. For example, we may expect a power law dependence close to r 2 / 3 for the 
structure functions. Direct numerical simulations have never shown a scale dependence for 
the dilatation and we may expect that it behaves like a relatively small factor. Then S(r) 
will still modify e as explained in the discussion above, however the power law dependence 
in r would be now slightly different. In conclusion and according to this simple analysis we 
see that compression effects (through the dilatation) will mainly impact the scaling law at 
the largest scales. 

Compressible spectrum. We may try to predict a power law spectrum for compressible 
turbulence. First, we note that several predictions have been made for the kinetic energy 
spectrum and also for the spectra associated with the solenoidal or the compressible part of 
the velocity 21|. We recall that although these decompositions are convenient for analytical 
developments, the associated energies are not inviscid invariants and the predictions are 
heuristic. For incompressible turbulence the situation is different because a prediction in 
f or ki ne tic energy spectrum may be proposed by applying a dimensional analysis 
directly on the 4/5s law [3j. Although it is not an exact prediction, the 4/5s law gives a 
stronger foundation to the energy spectrum for which a constant flux is expected. This re- 
mark was already noted in particular in recent 3D direct numerical simulations of isothermal 
turbulence where it is observed that the Kolmogorov scaling is not preserved for the spectra 
based only on the velocity fluctuations [15|. 

We shall derive a power law spectrum for compressible turbulence by applying a dimen- 
sional analysis on equation (ITT1) . Dimensionally, we may find e e gr ~ pu 3 . By introduc- 
ing the density-weighted fluid velocity, v and following Kolmogorov we obtain 
E v (k) ~ elgk~ 5 / 3 , where E v (k) is the spectrum associated to the variable v. Our predic- 
tion is compatible with the measurements recently made by direct numerical simulations 
[ljj] where the authors have noted that the exponent of the third-order velocity structure 



S 



mg 



2l| 



function is close to one if the field used is v instead of u. (Note that two other sea 
relations may be predicted like for the pressure term.) As explained by several authors 
in compressible turbulence we do not expect a constant flux in the inertial range. Here, 
the same conclusion is reached since we are dealing with an effective mean energy transfer 
rate. More precisely if we expect a power law dependence in k for the effective transfer 
rate one arrives at the conclusion that a steeper power law spectrum may happen at the 
largest scales. According to relation ([16]) and the simple estimate, pv 2 ~ r 2 / 3 ; we could 
have E v (k) ~ k~ 19 / 9 . This prediction means that for a small prefactor in ([16]) one needs an 
extended inertial range to feel the compressible effects on the power spectrum. The scale 
at which the transition happens between —19/9 and —5/3 may be the sonic scale k s as 



proposed in 



23] where such power laws were detected; in our case, a rough estimate gives 



fc s ~((V-u)/6u). 

Conclusion. The present work opens important perspectives to further understand the 
nature of compressible turbulence in the asymptotic limit of large Reynolds numbers with 
the possibility to extend the analysis to magnetized fluids with possibly other types of 
closures (e.g. polytropic gas), or to improve intermittency models by using the new relation 
- obtained by a statistical analysis at low order - as pivotal for a heuristic extension to 
statistical laws at higher order. We believe that astrophysics (e.g. interstellar turbulence) is 
one of the most important domain of application of the present work flfj| . 
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